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9P ■ Abstract 

G^ 

I This work presents a newly renovated approach to the analysis of second-order 

Riccati equations from the point of view of the theory of Lie systems. We show 
that these equations can be mapped into Lie systems through certain Legendre 
transforms. This result allows us to construct new superposition rules for studying 
second-order Riccati equations and to reduce their integration to solving (first-order) 
Riccati equations. 

^ ■ Keywords: Lie systems; superposition rules; second-order Riccati equations; Riccati 



equations. 

1 Introduction 

A big deal of work has lately been devoted to the analysis of the so-called Lie systems 
[2]-[ll]. Such systems are characterized by admitting superposition rules, i.e. functions 
which determine their general solutions in terms of generic families of particular solutions 
and a set of constants. Lie systems also enjoy a plethora of geometric properties [5]. This 
makes them suitable for the investigation of many systems of differential equations and 
relevant problems in the physics, mathematics and control theory (see [H] and references 
within) . 
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The aim of this work is to present a new Lie systems approach to second-order Ric- 
cati equations. Such equations, whose interest rehes on their appearance in physics and 
mathematics hterature [12]- [15], had aheady been investigated from the viewpoint of Lie 
systems in several works [6]- [8]. All these papers were based on the use of ad hoc changes 
of variables as well as some recent geometric techniques. However, we here introduce 
a further geometric method for the study of such equations, which is simpler and more 
powerful than previous ones. As a result, we derive a new superposition rule for analyzing 
second-order Riccati equations and prove that their explicit integration reduces to solving 
(first-order) Riccati equations [TT] . 



2 Fundamentals on Lie systems 

For simplicity, we generally assume all functions to be real, smooth, and globally defined. 
In addition, we restrict ourselves to analyzing problems on linear spaces. In spite of 
this, our results can easily be generalized and proved in full detail through several minor 
technical modifications. 

A t-dependent vector field over M" is a map X : (t, x) G M x M" t-)- X{t, x) G TW^ such 
that ToX = 112, where tts : {t,x) eRxW x e M" and r : TM" M" is the projection 
related to the tangent bundle TM". 

We call integral curves of a time- dependent vector field X over integral curves of 
the distribution spanned by its suspension, i.e. the vector field X = d/dt + X{t,x) over 
M X [1]. Given an integral curve 7 : s G M H- (t{s),x{s)) G M x M" for X passing 
through the point (to^^^o) G M x M", there exists a time-reparametrization i = i{s) such 
that 

which is referred to as the associated system related to X. Note that every t-dependent 
vector field determines a system of first-order differential equations of the above type 
and vice versa. This justifies to use X to denote both a t-dependent vector field and its 
associated system. 

Definition 1. A superposition rule for a system of first-order differential equations 

d/X^ 

— =X\t,x), ^ = l,...,n, (1) 

is a function $ : M"™ x M" — > R", of the form x = . . . , ki, . . . , kn), such that 

the general solution x(t) of ([T]) can be written as x{t) = $(x(i)(t), . . . , a;(m)(t); ki, . . . , kn), 
where X(i)(t), . . . ,X(m)(t) is any generic family of particular solutions and ki . . . ,kn is a 
family of constants to be related to initial conditions. 

Systems admitting a superposition rule are called Lie systems. The conditions guar- 
anteeing that the system ([1]) admits a superposition rule are stated by the Lie-Schejfers 
Theorem [2], Theorem 44]. A modern statement of this relevant result is described next 
(we refer to [U Theorem 1] for details). 
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Theorem 2. (Lie— Scheffers Theorem) A system (CP admits a superposition rule if 
and only if the associated t- dependent vector field X can be brought into the form Xt = 
Yla=i ba(t)Xa, for a certain family bi(t), . . . , Kit) of t-dependent functions and a set of 
vector fields Xi, . . . , Xr spanning an r-dimensional Lie algebra, the so-called Vessiot- 
Guldberg Lie algebra. 

The Lie-Scheffers Theorem enables us to reduce the integration of a Lie system to 
solving a special type of Lie system on a Lie group. In fact, every Lie system X on 
M" associated with a Vessiot-Guldberg Lie algebra V, e.g. Xt = J2a=i^a(t)Xa with 
Xi . . . ,Xr a basis for V, can be associated with a (generally local) Lie group action 
$ : G X — )• whose fundamental vector fields coincide with those of V. This action 
allows us to bring the general solution x{t) of X into the form x{t) = ^{g{t),Xo), where 
Xq G M" and g{t) is the solution of the Lie system 



with g{0) = e and X^, . . . , X^ being a basis of right-invariant vector fields over G closing 
on opposite commutation relations to Xi, . . . , X,, (see [S]). Conversely, the general solution 
of X determines the solution for Eq. ([2]) [HI [ID] • 

Several methods can now be applied to solve Eq. ([2]) and, consequently, the Lie system 
([1]). If g ^ TgG is solvable, the theory of reduction of Lie systems devised in [10] allows us 
to solve ([2]) by quadratures. More generally, we can use a Levi decomposition of g to write 
g ~ r©s (si © ■ ■ ■ © S;), where r is the radical of g, 5i © ■ ■ ■ ©S/ is the direct sum of a family 
of simple Lie subalgebras of g, and ©^ denotes a semi-direct sum of r and © ■ ■ ■ ©S/. The 
theorem of reduction of Lie systems [TOl Theorem 2] yields, with the aid of the previous 
decomposition, that the solution of a Lie system like ([2]) but defined on a Lie group with 
Lie algebra Si © ■ ■ ■ © enables us to construct a t-dependent change of variables that 
transforms Eq. ([2]) into a Lie system on a Lie group with Lie algebra r, which is integrable 
by quadratures (see [TU] for details). Summarizing, the explicit integration of ([2]) reduces 
to providing a particular solution of a Lie system related to Si © ■ ■ ■ © S;. 

General solutions of Lie systems can also be investigated through superposition rules. 
There exist different procedures to derive them [31 HI [9], but we hereafter use the method 
devised in [1], which is based on the diagonal prolongation notion [H [TT]. 

Definition 3. Given a t-dependent vector field X over W^, its diagonal prolongation X 
to M"(™+^) is the unique t-dependent vector field over M'^('^+i) such that: 

• Given pr : (x(o), . . . , X(^)) G M-(-+i) ^ X(o) G M", pi.Xt = XtWte M. 

• X is invariant under the permutations O with i,j = 0, . . . , m. 

The procedure to determine superposition rules described in ^ goes as follows. Take 
a basis Xi, . . . , Xr of a Vessiot-Guldberg Lie algebra V associated with the Lie system. 
Choose the minimum integer m so that the diagonal prolongations to M"™ of the elements 




(2) 
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of the previous basis are linearly independent at a generic point. Obtain n functionally in- 
dependent first-integrals Fi, . . . , F„ common to all the diagonal prolongations, Xi, . . . , Xr, 
to M'^('"+^), for instance, by the method of characteristics. Assume that these integrals take 
certain constant values, i.e., Fi = ki with i = 1, . . . ,n. Employ these equalities to express 
the variables of one of the copies of M" within ^"("^+1) in terms of the remaining variables 
and the constants ki, . . . ,kn- The obtained expressions constitute a superposition rule in 
terms of any generic family of m particular solutions and n constants. 

3 A new Lie systems approach to second-order Ric- 
cati equations 

The most general class of second-order Riccati equations is given by the family of second- 
order differential equations of the form 

+ ifoit) + fiit)x)— + Co(t) + ci(t)x + C2it)x^ + C3(t)x3 = 0, (3) 



with 



These equations arise by reducing third-order linear differential equations through a di- 
lation symmetry and a time-reparametrization [12]. Their interest is due to their use in 
the study of several physical and mathematical problems P, [3, El [El [I3] . 

It was recently discovered that every second-order Riccati equation ([3]) admits a t- 
dependent non-natural regular Lagrangian of the form 

L{t, X, v) = 



V + f/(t, x) ' 

mth U(t,x) = ao{t) + ai(t)x + a2(t)x'^ and ao(t) , ai{t) , a2{t) being certain functions related 
to the t-dependent coefficients of ([3]), see [12]. Therefore, 



^~ dv {v + u{t,x)y' ^ ' 

and the image of the Legendre transform FL : {t,x,v) G W C M x TM i— )■ {t,x,p) G 
M X T*M, where W = {{t,x,v) e R x TR \ v + U{t,x) ^ 0}, is the open submanifold 
M X where O = {{x,p) G T*]R | p < 0}. The Legendre transform is not injective, 
as {t,x,p) = ¥L(t,x,v) for v = ±1/ — U{t,x). Nevertheless, it can become so by 
restricting it to the open W+ = {(t,x,v) G M x TM | v + U(t,x) > 0}. In such a case, 
V = 1 1 ^—p — U(t,x) and we can define over M x O the t-dependent Hamiltonian 

h{t, x,p) = p ( — U{t,x)] — = —2y/^ — pU{t,x). 



'-P 
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Its Hamilton equations read 



U{t,x) = — ao(t) — ai{t)x — a2{t)x 



dt dp (^5) 
dp dh dU , , / / x „ / X X 

Since the general solution x{t) of every second-order Riccati equation (jlj) can be recov- 
ered from the general solution (x(t),p(t)) of its corresponding system ([5]), the analysis of 
the latter provides information about general solutions of second-order Riccati equations. 

The important point now is that system (|5]) is a Lie system. Indeed, consider the 
vector fields over O of the form 

V V ^ V ^ ^ V 2 n ^ 

Ai = — , A2 = A 3 = X- p—, X4 = X 2xp— , 

y/—p OX OX OX op OX op 

X d ^ / — d 

a/— p OX Op 

Their non- vanishing commutation relations read 



[XijXs] — 2^1' [-^15-^4] — ^5; [-^25^3] — ^2, [-^25-^4] — 

[-^2,^5] = -^1, [-^35-^4] = -^4, [X35X5] = 2^^' 



(6) 



and therefore span a five-dimensional Lie algebra V of vector fields. Additionally, the 
t-dependent vector field Xt associated with (|5]) holds 

Xt=Xi- ao{t)X2 - ai{t)Xs - a2{t)X^. (7) 

In view of expressions ([6]) and (j7]), system ([5]) is a Lie system. Note also that a similar 
result would have been obtained by restricting the Legendre transform over the open 
W_ = {(t, X, e M X TM I w + U{t, x) < 0}. 

Let us use the theory of reduction of Lie systems to reduce the integration of to 
solving a Lie system on a Lie group. Using a Levi decomposition of V, we get V ^ 
Vi (Bs V2, with V2 = {X2, X3, X4) being a semisimple Lie algebra isomorphic to sI(2,M) 
and Vi = (XijX^) the radical of V. Hence, V is isomorphic to the Lie algebra of a Lie 
group G = X SL{2,M.), where x denotes a semidirect product of by S'L(2,M), and 
there exists a local action $ : G x (9 — )■ (9 whose fundamental vector fields are those of 
V. It is a long, but simple, computation to show that 



where x = {ax + /3)/ (7X + 6), p = p (7X -|- 5)^ and a6 — P'-f = 1, is one of such actions (for 
a detailed example of how to derive these actions see [TI] Ch. 2]). 
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The above action enables us to write the general solution ^{t) of system ([5]) in the 
form C,{t) = ^{g{t),^o), where C,o & O and g(t) is the solution of the equation 

^ = - {X,^{g) - ao{t)X^{g) - a,{t)X^ig) - a,{t)X^{g)) , g{0) = e, (8) 

on G, with the X^ being a family of right-invariant vector fields over G whose vectors 
X^{e) G TeG close on the same commutation relations as the X^ (cf. [TT]). 

We now turn to apply to Lie systems (|8]) the theory of reduction for Lie systems. 
Since TgG ~ (Bs s[(2,]R), a particular solution of a Lie system of the form ([8]) but 
over S'L(2,M), which amounts to integrating (first-order) Riccati equations (cf. [Sj fTT]). 
provides us with a transformation which maps system (|H]) into an easily integrable Lie 
system over M?. In short, the explicit determination of the general solution of a second- 
order Riccati equation reduces to solving Riccati equations. 

Another way of analyzing the solutions of ([5]) is based on the determination of a 
superposition rule. According to the method sketched in Sec. 2, a superposition rule 
for a Lie system which admits a decomposition of the form ([7]), can be obtained 
through two common functionally independent first-integrals for the diagonal prolonga- 
tions Xi, X2, X3, X4, X5 to a certain T*M("'+^) provided their prolongations to T*M''" are 
linearly independent at a generic point. In our case, it can be easily verified that m = 4. 
The resulting first-integrals, derived through a long but easy calculation (see |8] for a 
similar procedure), read 

^0 = (a;(2) -X(3))y/P[2]Pi3) + (a;(3) -a;(i))v/P(3)P(i) + (a;(i) - x^2))^Pii)Pi2), 

Fi = (x(i) -X(2))VP(i)P{2) + (x(2) -a;(o))VP(2)P(o) + (a;(o) - a;(i))v/P(o)P(i), 

^2 = (a;(i) -X(^3))^/mm + (3^(3) -a;(o))v^(3)P(oy+ (3^(0) - X{i))VPio)P{i)- 

Note that given a family of solutions (x(j)(t),p(j)(t)), with i = 0, . . . , 3, of ([5]), then 

dFj/dt = XtFj = for J = 0, 1,2 and F, = F,(x(o)(t),P(o)(t), • • . , X(3)(t),P(3)(t)). 

In order to derive a superposition rule, we just need to obtain the value of p(o) from 
the equation ki = Fi, where ki is a real constant, to get 

r- _ + (a^(2) -X(l))^P(l)P(2) 

^ ^^^^ {X{2) - X{0))y/-P(2) + {X(0) - X{l))y/-P(l)' 

and then plug this value into the equation k2 = F2 to have 

^ ^ feir(X(i),P(i),X(3),P(3)) + /C2r(x(2),P(2),a:(l),P(l)) - Fo3;(i)^-p(i) 

ki{v-P{i) - V-Pm) + k2{y/-P{2) - V-Pw) - V-Pw^o 

P{0) = - [A;i/Fo(^-p(3) - >y-P{i)) + A;2/i^o(V-P(i) - \/-P(2)) + V-P(l)J^ 

where r(x(i),p(j), X(j),p(j)) = ^— — \/—P{j)X{j). The above expressions give us a 
superposition rule $ : (a;(i),p(i), X(2),P(2), a:(3),P(3); A;2) G T*M^xM2 ^ (x(o),p(o)) G T*M 
for system (|5]). Finally, since every X{i){t) is a particular solution for (|3]), the map T = ro$ 
furnishes the general solution of second-order Riccati equations in terms of three generic 
particular solutions X(i)(t), X(2)(t), X(3)(t) of ([3]), the corresponding p(i)(t),p(2)(t),P(3)(t) 
and two real constants ki,k2. 
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4 Conclusions 



Our work introduces many advantages with respect to previous methods for studying 
second-order Riccati equations [n]-[H]- Apart from avoiding the use of ad hoc changes 
of variables as in [HI [7], we directly transform second-order Riccati equations into Lie 
systems (|5]) by Legendre transforms. This is much simpler than other approaches based 
on more elaborated geometric theories [8] , which map second-order Riccati equations into 
Lie systems associated with a 5 [(3, M) Lie algebra. Indeed, our method reduces the explicit 
integration of second-order Riccati equations to solving Lie systems related to s[(2,]R), 
e.g. Riccati equations. 
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